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Empirical evidence in heavy fermion, pnictide, and other systems suggests that un¬ 
conventional superconductivity appears associated to some form of real-space electronic 
order. For the cuprates, despite several proposals, the emergence of order in the phase 
diagram between the commensurate antiferromagnetic state and the superconducting 
state is not well understood. Here we show that in this regime doped holes assemble in 
“electronic polymers”. Within a Monte Carlo study we And, that in clean systems by 
lowering the temperature the polymer melt condenses first in a smectic state and then 
in a Wigner crystal both with the addition of inversion symmetry breaking. Disorder 
blurs the positional order leaving a robust inversion symmetry breaking and a nematic 
order, accompanied by vector chiral spin order and with the persistence of a thermo¬ 
dynamic transition. Such electronic phases, whose properties are reminiscent of soft 
matter physics, produce charge and spin responses in good accord with experiments. 

PACS numbers: 74.25.Ha, 71.28.+d, 75.25.-j 


Introduction 

The anomalous behavior of several physical quantities 
above the superconducting transition temperature sug¬ 
gests that high-temperature superconductivity emerges 
close to a quantum critical point between a broken- 
symmetry state and the disordered phase, in analogy 
with heavy fermion materials, pnictides and organics [T] . 
However, identifying the broken-symmetry phases has 
proven much more difficult than in other materials, de¬ 
spite several “gold rushes” in the underdoped region of 
the phase diagram, triggered by the observation of stripe 
order nia, nematic order w, time-reversal symmetry 
breaking [7], and incommensurate charge-density-wave 
order [5]. 

While some form of charge order (GO) is well es¬ 
tablished [9], an important question is how this order 
is formed starting from the two extremes of the phase 
diagram. When coming from the high-doping region, 
GO seemingly arises as a second-order instability of 
the uniform strongly correlated metallic state, produc¬ 
ing incommensurate charge density waves driven by mag¬ 
netic mm. phononic mm. or mixed |14j microscopic 
mechanisms. On the other hand, the occurrence of CO 
from the Mott insulating low-doping side is more directly 
tied to the tendency of Mott antiferromagnets to expel 
and segregate charges [T5HT8] .This latter region is the 
playground of our present work. 

Recently, it was pointed out [iniiiD] that, at very low 
doping, charge segregation may acquire features, related 
to the occurrence of topological excitations in doped anti¬ 
ferromagnets. The starting point is the observation that 
holes in an antiferromagnet induce a vortex (V) or an¬ 
tivortex (A) texture in the surrounding spin ordering. 
While isolated vortices are energetically expensive, a VA 


pair is stable EIHMI, because its annihilation is hindered 
by the strongly correlated character of the doped holes 
and the disturbance of the antiferromagnetic background 
rapidly dies out at large distances. As in early propos¬ 
als [25l |26] inspired by the work of Villain [27l [28] , these 
VA “dimers” explain the extremely rapid destruction of 
long-range antiferromagnetic order with doping. 

In the present scenario the dimers or “nematogens” 
self-organize and give rise to “electronic soft matter” ef¬ 
fects |29|. Specifically, the dimers may undergo a “poly¬ 
merization process”, triggering charge segregation into 
segments, tightly bound to V and A spin textures. These 
segments not only align forming a nematic state, but can 
also break inversion symmetry mm. due to their in¬ 
trinsic topologic dipolar character (associated with the V 
and A at the endpoints of the “polymer”). This state, 
which was named ferronematic miiDi, is accompanied 
by a spin spiral state sustaining a net spin current. At 
large scales, this feature is reminiscent of other propos¬ 
als (SUMSSj) which are however based on impurity states 
instead of the polymer states which are central to our re¬ 
sults. A ferronematic phase was proposed also to occur 
in ultra-cold dipolar Fermi gases of atoms |M| . 

We pose here the following fundamental questions: 
which other phases can be sustained by the electronic 
polymers, how are they affected by quenched disorder, 
what is the fate of the thermodynamic phase transitions 
expected in ideally clean systems, and how their char¬ 
acteristic temperature scales emerge from the (usually 
much higher) electronic scales of the system. In order to 
study the problem at the large length scales probed by ex¬ 
periments we carry out a multiscaling approach starting 
from a microscopic model and derive a mesoscale effec¬ 
tive model treated with Monte Carlo. We obtain a rich 
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FIG. 1. Charge and spin configurations in the dif¬ 
ferent phases. White and black circles represent the posi¬ 
tive and negative topological charges, respectively. The dif¬ 
ferent colours denote the angle of the staggered magnetiza¬ 
tion. The images are Monte Carlo snapshots in the absence 
of quenched disorder (a-c) in the thermally disordered phase 
with T = 50K (a), in the ferrosmectic phase at T = 38K (b), 
in the ferrocrystal phase at T = 8 K (c) and in the ferrone¬ 
matic phase at T = 40 K (d) which appears in the presence of 
quenched disorder (Qion/Qrep = 0.125) . The white lines in 
(b) highlight the “triangular” arrangement of the segments. 

phase diagram for the electronic polymers as a function 
of temperature and disorder which allows to rationalize 
the charge and spin responses observed experimentally. 


Results 

Numerical Simulations. We start from the very low 
doped regime of few holes in the spin antiferromagnetic 
background of a Cu02 plane modeled by a one-band Hub¬ 
bard model. We study the dimers at mean-field level in 
the Gutzwiller approximation (GA) (see Methods, Sup¬ 
plementary Note 1 and Supplementary Fig. 1). With 
realistic parameters for La 2 -xSra;Cu 04 , the most favor¬ 
able configuration for two holes is along the diagonal of a 
plaquette with a planar dipolar distortion of the antifer¬ 
romagnetic background. The latter can be visualized as 
due to a V and an A centered close to (but not exactly 
at) the vertices of the plaquette and forming a “topolog¬ 
ical dipole” (TD). There is another two-hole mean-field 
solution which is non-planar and consist of a skyrmion 
texture m which, for the present parameters, is ~ 100 K 
higher in energy than the TD and therefore will be ne¬ 
glected at low temperatures in the following. 


Studying metastable planar configurations, in which 
two or more of these TDs are arranged with different po¬ 
sitions and orientations (Supplementary Figs. 2 and 3) 
we find, as expected, that at large distances holes inter¬ 
act through a logarithmic interaction |35j between their 
topological charges, while at short distances their interac¬ 
tion is modified by quantum effects related to the overlap 
of the hole wave functions. The logarithmic interaction 
stems from the fact that for planar textures the long- 
range behavior can be captured by an XY model [551 US] ■ 
Notice that we are not claiming that the symmetry of the 
model is reduced from Heisenberg to XY. Indeed each 
texture has a zero mode related to the change of the 
plane which contains the spins, as it should for an 0(3) 
symmetric model. However, contrary to what would hap¬ 
pen for a single vortex in the pure Heisenberg model, the 
textures have no unstable modes [50] which would break 
the planar character of the texture, i.e. they are locally 
stable and their energy is correctly captured by a planar 
magnetic model. 

In order to enable simulations in large systems we do 
not consider explicitly the spin degrees of freedom but in¬ 
tegrate them out to generate effective interactions among 
topological charges. While this is an enormous compu¬ 
tational advantage, it limits our simulations to low dop¬ 
ing (rifi O.I holes per unit cell) where spin currents 
are small on average and the superposition principle is 
valid, allowing for a mapping of topological charges onto 
a two-dimensional (2D) Coulomb gas |35|. The effec¬ 
tive interaction among topological defects, needed for the 
Coulomb gas model, is obtained by fitting the energy 
of several metastable zero-temperature GA solutions ob¬ 
tained in the Hubbard model. 

As a consequence of the interaction mediated by the 
antiferromagnetic background, when a large even num¬ 
ber of holes is added to the system, these tend to bind 
into a single polymeric chain of alternating topological 
charges, ending with a V and an A. Adding the real 
three-dimensional long-range Coulomb repulsion among 
holes, whose strength is measured by a parameter Qrep 
(see Methods), these long polymers break into smaller 
polymers, as shown in Fig. [^ and Supplementary Fig. 5. 

We work at temperatures T smaller than the binding 
energy of individual VA pairs (« 100 K), so that the num¬ 
ber of unbound topological charges is negligible. There¬ 
fore, our basic constituents in the Monte Carlo compu¬ 
tations are the TDs. These are modeled by a bound V 
and A, each moving on the sites of a square lattice, with 
the topological charge adjusted so that the dipole mo¬ 
ment matches the Gutzwiller computations (see Meth¬ 
ods). Since there are no topological constraints on the 
charge ±fc of the V and the A, they turn out to be frac¬ 
tional, k « 0.8 (see Supplementary Note 1). 

A crucial problem in cuprates is to determine how dis¬ 
order affects the ordered phases of the ideal “clean” sys¬ 
tem [351 . In order to address this issue, the holes at- 
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FIG. 2. Charge and spin structure factor. Density plots 
in the 2D reciprocal space for: the charge structure factor for 
the clean system at (a) 38 K and (c) 8 K; the charge structure 
factor for the system with Qion/Qrep = 0.125 at (b) 38 K and 
(d) 8 K. The red solid circles represent the position of the spin 
peaks in the reciprocal space [shifted by qAp]- The arrow in 
panel (c) shows a ferrocrystal peak. The lower panels show the 
diagonal cut of the charge and staggered spin structure factor 
in the 2D reciprocal space at 38 K for (e) the clean system and 
(f) the system with Qion/Qrep = 0.125. In order to see more 
clearly the effects of the broken C 4 symmetry the averages are 
restricted to configurations with 4 > > 0 corresponding to the 
expected response in a single domain sample. 


tached to the topological charges are subject to a ionic 
disorder potential with strength Qiom generated by the 
counterions out of the CUO 2 plane (see Methods and 
Supplementary Note 2 and Supplementary Fig. 4). The 
magnitude of Qion is difficult to estimate because it de¬ 
pends on screening processes not comprised in the model. 
Therefore, we treat Qion/Qrep as a phenomenological di¬ 
mensionless parameter which characterizes the amount 
of disorder. 

We consider a L x L square cluster with even number 
of holes Nh, corresponding to Nh/2 TDs. Although we 
explored various Filings, for the sake of definiteness in 


the present Communication we report the results for the 
typical case L = 100 and = 300, corresponding to a 
hole doping Uh = 0.03. 

To characterize the broken symmetries, we define a 
nematic order parameter (/>(T') [see Eq. ([^ in Methods], 
which becomes different from zero when the C 4 rotational 
symmetry of the lattice is broken. We also define the 
polarization of the system as the normalized sum of all 
the TD moments projected on the (1,1) and (1,-1) pre¬ 
ferred directions [cf. Methods, Eqs. @ 0 )]. A nonzero 
polarization in the system implies a breaking of inversion 
symmetry of the magnetic texture. Vector chiral spin or¬ 
der is characterized by the chirality Xi,±i [see Eq. (|^ in 
Methods]. Einally, the charge and spin structure factors 
allow us to further characterize the various phases. 

The Monte Carlo computations find at high temper¬ 
ature a classical liquid of dimers, which tend to form 
longer polymers as temperature is lowered (Supplemen¬ 
tary Eigs. 6), and to align along the diagonal directions, 
which are energetically favorable. Eig. & reports a snap¬ 
shot of this high-temperature phase, taken during the 
Monte Carlo evolution. 


Clean System. For the clean system (Qion = 0) we find 
that, when T is low enough, the segments orient to form a 
state with C 4 symmetry breaking. As is clearly visible in 
Fig. (see also Supplementary Fig. 5), associating the 
segments with “polymers”, the low-T phase corresponds 
to the so-called smectic order of soft matter [37] , in which 
the system has long-range positional order in one direc¬ 
tion, with periodicity £c = l/^c [<lc being the magni¬ 
tude of a characteristic wave-vector in reciprocal lattice 
units (rlu)], but remains “liquid” in the other direction. 
This manifests as sharp (resolution limited) peaks in the 
structure factor along the diagonal of the Brillouin zone 
(Fig. ^), i.e., perpendicular to the preferred polymer 
direction, signaling long-range positional order. As is 
shown in panel (e), the main spin peak is twenty times 
higher than the charge peak, because the spectral weight 
of the latter is spread over a wider range of wave-vectors, 
due to the highly anharmonic charge distribution. (No¬ 
tice that both structure factors, as defined in Methods, 
have the same normalization.) Regarding CO, this state 
has the same symmetry as a diagonal stripe state. How¬ 
ever, the charge is uniform along the stripe direction only 
after thermal fluctuations have been taken into account. 
In addition, this state breaks inversion symmetry, i.e., 
TDs tend to point in the same direction, thus we call it 
ferrosmectic. 

The ferro ordering associated with this and other 
phases is not trivial. Indeed, in contrast to dipoles on 
a cubic lattice in three dimensions, it would not occur if, 
for example, the TDs were arranged at fixed positions on 
a square lattice. This stems from the 2D dipole-dipole 
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interaction, which is ferroelectric in a nearly head-to-tail 
configuration of the dipoles, but is antiferroelectric for a 
side-by-side configuration. In our model, the ferro ten¬ 
dency wins because the real Coulomb interaction between 
the electrically charged holes favors short-range triangu¬ 
lar arrangements of the segments, i.e., segments in one 
row tend to face gaps in the neighbouring rows as is clear 
in Fig. (highlighted by white lines) and c, so that the 
side-by-side arrangements are rare. The colours in Fig. 
show the phase of the local staggered magnetization. In 
panels b-d the phase increases monotonically along one 
diagonal, indicating that these phases have long-range 
vector chiral order, i.e., Xi,-i 7 ^ 0 or xi,i 7^ 0. 

Upon further lowering the temperature, the ferrosmec- 
tic phase keeps the ferro ordering (and the vector spin 
chirality), but forms a Wigner crystal for T < lOK as 
shown in Fig. [^. This “ferrocrystal” manifests as addi¬ 
tional resolution limited off-diagonal peaks in the charge 
structure factor (as the one indicated with an arrow in 
Fig. [^) and which again signal long-range charge and 
spin order. 

Effect of Disorder. The properties of the phases 
change dramatically upon the introduction of quenched 
disorder. The ferrosmectic and ferrocrystal peaks 
broaden and weaken very rapidly (Fig. ,d,f), thus 
long-range positional order is lost and the ferrosmectic- 
ferrocrystal transition is smeared. Remarkably long- 
range nematic and vector chiral order (accompanied by 
inversion symmetry breaking) remain at finite disorder, 
and the phase becomes the ferronematic state proposed 
in Refs. [laiig. Fig. shows how CO is almost en¬ 
tirely destroyed by small disorder. Long range spin or¬ 
der also is destroyed but short-range spin order, sig¬ 
naled by incommensurate peaks, whose width is well re¬ 
solved in our system size, persists. For Qion/Qrep > 0.25 
even the broad incommensurate magnetic peaks disap¬ 
pear (Supplementary Note 4 and Supplementary Figs. 
6-8). This would contradict experiments, we thus esti¬ 
mate Qion/Qrep < 0.25 in real systems. 

The red solid circles in Fig. -d show the vectors 
±ie/V2,-e/V2)T.\ .u. where e is the magnetic incom¬ 
mensurability, i.e., for the orientation of Fig. [^,c mag¬ 
netic peaks appear at qaf i (e/-\/2, —e/v^) r.l.u., with 
QAF = (0.5,0.5) r.l.u. From all panels we see that the 
main magnetic peaks appear at half the incommensurate 
wave-vector of the main charge peaks. At first sight, this 
relation, well known for spin collinear stripes [38l [39] , is 
surprising here, since the incommensurability should be 
linked to the topological polarization [12] . However, close 
inspection of Fig. reveals that each segment acts as 
an antiphase domain wall for the antiferromagnetic back¬ 
ground, yielding jumps of the phase of the magnetic order 
parameter close to tt upon crossing the line of polymers. 
On the other hand, the phase is approximately constant 
in between two polymer rows. Thus, the magnetization 
behaves similarly to the case of a collinear stripe array. 



q[rlu] r[K] 

FIG. 3. Commensurate-incommensurate transition 

(a) Diagonal cuts of the spin structure factor for different 
temperatures as a function of momentum with q defined as 
in Fig. and disorder Qion/Qrep = 0.125. (f). The peaks 
have been convoluted with a Gaussian (standard deviation 
0.041 [rlu]) to take into account a finite experimental resolu¬ 
tion. (b) Height of the structure factor shown in (a) at the 
commensurate antiferromagnetic wave-vector (blue) and at 
the incommensurate position with respect to the background 
(red) as a function of temperature. The green data (right 
scale) shows the incommensurability as a function of tempera¬ 
ture. The vertical line marks the ferronematic transition. The 
inset shows the experimental peaks height from Ref. |40| for 
doping Uh — 0.0192, slightly below the complete disappear¬ 
ance of static antiferromagnetic order as reveled by muons. 
The evolution of the incommensurate peaks has been shown 
to be continuous m across the critical doping nu = 0.02. 


Spin canting produces small corrections to the ’factor of 
two’ relation, which are below our momentum resolution 
to be visible in Fig. 

Raising the temperature at small disorder, the broad¬ 
ened spin and charge peaks gradually decrease (Fig. i^) 
without any sign of a sharp transition in the intensity 
(Fig. I and Supplementary Fig. 6, respectively) as also 
observed experimentally at similar dopings (see inset of 
Fig.[3]D and Ref. m)- In contrast, studying the polar¬ 
ization and nematic order parameter distribution we find 
that the transition from the ferronematic to the melted 
polymers is of first order and remains sharp for our sys¬ 
tem size (see Supplementary Note 5). Thus a thermody¬ 
namic transition persists even in the presence of disor¬ 
der. The thermodynamic transition temperature is sig¬ 
naled by a change of behavior in the magnetic structure 
factor from commensurate to incommensurate, providing 
a simple experimental tool to detect the transition line 
(Fig. ^). This is because the incommensurability e is 
related to degree of polarization in the system and thus 
acts as an order parameter [12] . 

Phase Diagram. Fig. [^ reports the phase diagram 
obtained from the above analysis. The ferrocrystal (thick 
yellow line) and ferrosmectic (thick pink line) phases are 
well defined only in the absence of disorder. At finite dis¬ 
order, they survive as short-range ordered states. This is 
indicated by the yellow shade for the ferrocrystal and by 
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the magenta shade for the ferrosmectic state. The light 
blue region is the long-range ordered ferronematic state, 
while the red line indicates the first-order transition to 
a liquid of short polymers. We never find a purely ne¬ 
matic phase, characterized by a nonzero nematic order 
parameter but zero polarization and zero global vector 
chiral spin order. The last phase, which is allowed by 
our model, could possibly be stabilized in a different pa¬ 
rameter regime, as an intermediate phase between the 
ferronematic and the disordered phase. 

Our results are in good qualitative agreement with the 
phase diagram obtained by completely different methods 
in Ref. |36] . On the other hand, we find an additional 
inversion symmetry breaking and we provide realistic es¬ 
timates of the parameters of the model, of the experi¬ 
mentally measurable structure factors and of the charac¬ 
teristic temperatures of the transitions. 

The order of magnitude of the transition temperature 
to a polarized state can be estimated using a mean-field 
approximation. For dipoles in two-dimensions at random 
positions but with a non zero average dipole moment (p), 
the dipolar field can be computed using elementary elec¬ 
trostatics to be Ed = 2 'k'^ ps{p)nh/Nc where ps is the 
magnetic stiffness of the system and the number of 
charges per segment. Assuming a mean-field approxima¬ 
tion where the dipoles, of strength pq = y/2{Nc — l)fc, 
can fluctuate in 4 possible orientations (f, —>■, i, ■<—) we 
obtain the ordering temperature, 

ksTc = {npof ps^ ■ ( 1 ) 

With the present parameters and assuming A"c ~ 4 (Sup¬ 
plementary Fig. 6), we find « 196 K. This is of the 
correct order of magnitude (he much smaller than the 
original electronic scales) taking into account that we 
have neglected the positional entropy which will reduce 
Te. 

Both ferrosmectic and ferrocrystal charge orderings are 
not commensurate. Thus, they break a continuous [U(l)] 
symmetry in two dimensions. Even in the presence of in¬ 
finitesimal disorder strict long-range order is forbidden 
[42j and one finds quasi long- or short-range order. In 
contrast, the nematic order parameter breaks a discrete 
(Z 2 ) symmetry and is much more robust against disor¬ 
der. In our computations we have in addition vector chi¬ 
ral spin order (or equivalently a topological polarization) 
which also breaks a discrete (Z 2 ) symmetry, but which 
does not couple linearly to the local disorder, in con¬ 
trast to the nematic order parameter mm]. General 
arguments indicate that the discrete symmetry breaking 
should be much more robust than the breaking of a con¬ 
tinuous symmetry [36l|45]5 we indeed find. We expect 
the nematic order to behave similarly to the random field 
Ising model: lacking long-range order in a strictly 2D sys¬ 
tem, but ordered within a correlation length which can 
be exponentially large for small disorder [46) . favoring 
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FIG. 4. Phase diagram as a function of temperature and dis¬ 
order strength. The yellow (pink) thick line at zero disorder 
corresponds to ferrocrystal (ferrosmectic) long-range order. 
The yellow region is short-range ferrocrystal order while the 
magenta region corresponds to the short-range ferrosmectic 
order. At finite disorder, below the red line, the system has 
long-range ferronematic order (light blue region) while a poly¬ 
meric liquid is found above the red line, up to the highest 
temperatures reached in our study. 

a crossover to three-dimensional long-range order in the 
presence of a small inter-layer coupling mm- 

Since the ferronematic state has short-range spin order 
and long-range vector chiral order (at q = 0 wave-vector), 
it can be identified with the chiral spin liquid believed to 
take place in frustrated magnets [4^50] . 


Discussion and conclusions 

Our results allow us to rationalize several experimental 
findings, and imply some predictions which have not yet 
been tested. 

Experiments show that hole doping destroys commen¬ 
surate antiferromagnetic order much more rapidly than 
what would be expected by site dilution mm- Fig.[ 2 f 
shows that this is explained by a small density of TDs. 
The ability of VA pairs to rapidly depress commensu¬ 
rate ordering was noticed before [IS] , although these 
authors did not consider the collective ordering of the 
dipoles. 

Incommensurate spin scattering has been detected in 
the early days of high-T^ [SS] and interpreted in terms of 
stripes [11[3|. However, stripes are associated with charge 
modulations which are extremely hard to measure, in 
contrast to spin modulations. CO generally emerges as¬ 
sociated with a structural distortion close to Uh = 1/8 
which can be controlled by codoping with Nd Hi or 
doping/codoping with Ba [54H55] . All these observa¬ 
tion of CO are at doping close to rih = 1/8. The in- 
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tensity of CO decreases strongly with underdoping and 
extrapolates to zero around nt ~ 0.09 [56]. To the best 
of our knowledge, incommensurate static charge order 
has never been reported in the present, heavily under¬ 
doped, regime, in contrast to incommensurate spin or¬ 
der m, which persists. This dichotomy is explained 
by our simulations which, while reproducing the incom¬ 
mensurate spin ordering, show very weak charge-ordering 
peaks, barely emerging from the background noise, even 
for weak disorder (Fig. [^). 

Close to Uh = 1/8 magnetic Bragg peaks appear quite 
sharp and often resolution limited [a ng l58] indicat¬ 
ing long-range order. As doping is reduced static peaks 
are still observed but become broad with a well resolved 
width of the order of the incommensurability indicating 
a correlation length of the order of the spin periodic¬ 
ity |561157| . This is in excellent agreement with our mag¬ 
netic structure factor in Fig. [^. We interpret this fea¬ 
ture as an indirect signature of long-range vector chiral 
spin order without long-range magnetic order i.e. the 
ferronematic state we propose. For rih = 0.03, experi¬ 
mental magnetic peaks have been detected with incom¬ 
mensurability e « 0.032 [41j in good agreement with our 
computations yielding e « 0.028 at low temperature. 

Neutron scattering experiments in Y-based materials 
have shown HE] that the magnetic incommensurability 
as a function of temperature behaves as an order pa¬ 
rameter. Such a behavior is naturally explained by our 
model, where the incommensurability, in the presence of 
weak or no CO, is closely linked to the topological po¬ 
larization dS], which is an order parameter (see Fig. |^. 
Furthermore we propose that the temperature at which 
the static magnetic structure factor changes from a dou¬ 
ble peak structure to a single peak structure is a proxy 
of the thermodynamic critical temperature below which 
long-range chiral spin order is established. 

The transition from incommensurate behavior to com¬ 
mensurate behavior has been observed also in the spe¬ 
cific La-family we focus on in the present computa¬ 
tions. Indeed experimental low-energy inelastic neutron 
scattering peaks as a function of temperature reported 
in Fig. 5a of Ref. [59] (see also Ref. (S^) show the 
same behavior as we find for the static structure fac¬ 
tor. However the transition from two incommensurate 
peaks to an antiferromagnetic commensurate peak takes 
place around 55 — 100 K. On the other hand quasistatic 
scattering shows a transition at around 20 — 30 K (see 
Refs. |40| im ]60] and inset of Fig. Our computa¬ 
tions provide an energy integrated structure factor which 
is expected to show the transition between the inelastic 
and quasistatic cases. Indeed we find the commensurate- 
incommensurate transition at around 45 K fully consis¬ 
tent with the neutron scattering measurements. Such 
an agreement on the temperature scales and qualitative 
behavior further supports our identification of the low 
temperature state observed in cuprates as a long-range- 


ordered ferronematic. 

Notice that, in contrast with the small ordering scales 
we find, the starting point electronic Hamiltonian has 
bare electronic scales of the order of 3000 K or more. This 
strong reduction of energy scales indicates that our multi¬ 
scale modeling has identified the correct dynamical vari¬ 
ables of the problem. Eq. Q shows that the energy scale 
is set by the magnetic stiffness and the density. Notice 
also that the proposed thermodynamic transition occurs 
at a temperature much lower than the pseudogap tem¬ 
perature (« 300 K) which instead nearly extrapolates to 
the Neel temperature of the undoped sample m- 

At even lower temperatures of the proposed ferrone¬ 
matic transition a so called cluster spin glass state is 
observed consisting of strongly coupled clusters of spins 
with weaker coupling among clusters EIlISlISnilSMBl]- 
The ferronematic state of Fig. corresponds precisely 
to this physical picture. 

Eq. Q predicts a linear relation between doping and 
the ordering temperature. However for nh > 0.02 one 
should take into account that the finite magnetic corre¬ 
lation length, the expected weakening of the stiffness by 
doping and the breakdown of the linear regime for the 
topological charges will lead to a slowing down of the dop¬ 
ing dependence. Interestingly for hole content nu < 0.02 
the temperature at which the cluster spin glass state is 
observed is linear with doping |5^HB4] consistent with our 
proposal. Numerically we find that the temperature of 
the ferronematic transition increases approximately lin¬ 
ear with doping as Tc ^ 1500 K compared with the 
experimental behavior Tc ~ 815 K Uh- Our larger sloop 
may be due to an overestimation of Ps, a smaller Nc 
and/or dynamical effects which may give an apparent 
shift of the transition. 

In the presence of spin-orbit coupling, long-range vec¬ 
tor chiral spin order gives rise to a real electric polariza¬ 
tion, i.e., the system becomes an improper ferroelectric 
[55] . Unfortunately, this effect is hard to observe be¬ 
cause, as soon as the system becomes metallic, it cannot 
support a finite electric polarization. Notwithstanding, a 
finite ferroelectric polarization has been reported at low 
temperatures in oxygen [66| and Li m doped La 2 Cu 04 , 
the samples having a strongly insulating character. The 
fact that the effect appears independently of the dopant, 
and that the remnant polarization can be oriented along 
different axes with external fields, clearly points to a 
magnetic origin of the ferroelectric polarization. Further¬ 
more, more recent experiments show a clear correlation 
between magnetoelectric effects and stripe orientation in 
Sr doped La 2 Ni 04 , suggesting that stripe effects are in¬ 
volved [SS]. Experiments at finite frequencies suggest 
that inversion symmetry breaking sets in at tempera¬ 
tures higher than the temperatures at which the sam¬ 
ple is insulating enough to support a static polarization. 
All these experiments support our conclusion that under¬ 
doped cuprates show long-range vector chiral spin order. 
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A possible test to our model would require second har¬ 
monic generation to detect inversion symmetry breaking 
in non-insulating samples. We predict that in the fer- 
ronematic phase the inversion symmetry breaking should 
track the behavior of the incommensurability as a func¬ 
tion of temperature. This relation, however, will break 
down in the collinear stripe phase found around nu = 
1 / 8 . ^ 

With the present method we cannot access quantita¬ 
tively the crossover to collinear stripes. In this regime, 
the mapping to the Coulomb gas breaks down due to non¬ 
linear effects. However, one can anticipate that the aver¬ 
age length of the segments will keep growing with dop¬ 
ing, leading to a concomitant increase of the ferrosmectic 
correlation length. According to our findings, the disor¬ 
der induced by the dopants will partially counteract this 
increase, but the associated impurity potential will also 
be progressively screened, opening the possibility that 
segments coalesce into stripes with long-range order and 
narrow magnetic peaks. 

We thus propose that underdoped cuprates have a 
long-range broken symmetry state at low doping. This 
puts the cuprate phase diagram into the same class of 
phase diagrams of a wide class of materials [1] in which 
unconventional superconductivity emerges from a phase 
characterized by real-space electronic long-range order. 


Methods 

Model. Treating the single-band Hubbard model 
within a Gutzwiller approximation, a single hole in the 
antiferromagnetic background is found to form a spin po- 
laron, while two holes tend to occupy the cores of a spin 
V and A that attract each other, thereby lowering their 
energy. The long-range part of this texture is treated 
using generalized elasticity m and exploiting the cor¬ 
respondence between a spin vortex and a 2D Coulomb 
charge [35] . In the absence of disorder and holes the mag¬ 
netic correlation length is expected to be very large but 
finite due to thermal fluctuations. This provides a natu¬ 
ral cutoff at a distance A for the long-range interactions 
between topological charges at large distances [Ss] 1^ . 
Therefore, the interaction energy is well described by 

n27T poo Aqr COS 9 

Virijr) = pskik 2 / d6> / dg 

JO Jo g + A 

= 2TTpskik2KQ , ( 2 ) 

where fci .2 = are the topological charges (in our case, 
k = 0.8) and Kq is the zeroth order modified Bessel 
function, which reproduces the logarithmic interaction 
at short/intermediate distances (r < A) and decays ex¬ 
ponentially at long distances (r ;g> A). We have checked 
that changing A the results are substantially the same, 
as long as it remains larger than the typical distance be¬ 
tween the segments. In the simulations we take A of the 


order of the system size for numerical convergence pur¬ 
poses. 

While Eq. ([^ reproduces well the energy of Gutzwiller 
calculations for the single band Hubbard model at large 
distances, as expected, it fails at short distances where 
the short-range physics of the Hubbard model becomes 
relevant. Therefore, we also include short-range terms 
extrapolated from the Gutzwiller calculations (Supple¬ 
mentary Note 1). 

Furthermore, each topological charge arising from the 
spin texture corresponds to a positive electrically charged 
hole in the Cu02 planes of the doped cuprate. Therefore, 
our model includes also the three-dimensional Coulomb 
repulsion potential which, for two holes at a distance r, 
is parametrized as 

= (3) 

r 

Here Qrep, incorporating, e.g., the static dielectric con¬ 
stant, represents the strength of the repulsion and is an¬ 
other parameter of our model. We fix Qrep/a = 49meV, 
where a is the in-plane lattice constant, so that the aver¬ 
age number of holes in a polymer, for very low density is 
Nc « 2. As the density increases, this number tends to 
increase too m, yielding the results of Supplementary 
Fig. 6 for the present density {rih = 0.03). 

Finally, the charged holes doped into the Cu02 planes 
leave back negative countercharges. For instance, in 
La 2 -a;Sr 2 ;Cu 04 , which we take as a prototype cuprate, 
negative Sr ions randomly replace La atoms between 
two consecutive planes. We therefore introduce disorder, 
generating a random distribution of point-like negative 
charges, one for each positive hole in the plane, which 
act as pinning centers for the carriers in the plane (Sup¬ 
plementary Note 2). The ions are located out of plane, 
at a distance d « 0.58a, from the center of the in-plane 
plaquette. Each impurity interacts with the holes in the 
plane through an attractive three-dimensional Coulomb 
potential 


where d is the distance between the hole and the impu¬ 
rity, and the strength of the interaction Qion measures 
the intensity of disorder. This procedure thus produces 
a disordered potential in the plane in which holes and 
their associated topological charges move. We show one 
realization of the impurity potential in Supplementary 
Fig. 4. 

Characterization of the phases. To characterize 
C 4 rotation symmetry breaking, we introduce the ne¬ 
matic order parameter 

^ ^(n(r,)n(rj + x + y) 

- n{r, + y)n{Ti + x)), 


(5) 





where n(ri) is the number of holes on the site labelled 
by Ti, and x,y are unit vectors along the corresponding 
directions of our 2D square lattice. The angular brackets 
imply thermal average, denotes the total number of 
charges, and the sum runs over all the lattice sites. To 
characterize inversion symmetry breaking, we introduce 
the polarization P = (Pa;,Pj,), as the normalized sum of 
all the TDs. Since in our model diagonal polarizations 
are favoured, we introduce the components 

P(l,l) = (P. + Py)/V2, (6) 

F(i,-i) = (F.-Py)/V2. (7) 

Vector chiral spin order is characterized by the parameter 


Xi,±i 


1 


^([s(ri) X s(r, + i±y)] • z), 


( 8 ) 


where s(ri) is the local spin density. 

Our Monte Carlo calculations also yield the thermal 
averages of the static charge (c) and spin (s) structure 
factors, S'c(q) = (1 - (5q,o)|7Xc(q)P/T^ and 5's(q) = 
|Ks(q)p/L2^ where 


F^c(q) 

K.(q) 


;^J^exp(iq-ri) s(ri). 


(9) 

( 10 ) 


With these definitions and using unit-length spins, the 
structure factors have the same normalization and satisfy 
Eq>5's.c(q) 

The ferronematic-ferrosmectic crossover in Fig. was 
characterized analyzing the height of the main charge 
peak as a function of temperature. 

Monte Carlo analysis. We carried out Monte Carlo 
calculations exploiting the parallel tempering technique. 
To analyze the spin degrees of freedom for a given config¬ 
uration, we attach to each topological charge the struc¬ 
ture of a (anti)vortex in the spin background and we 
perform a linear superposition, allowing then each spin 
to relax according to the XY Hamiltonian [35]. Sup¬ 
plementary Fig. lb reports an example for the case of 
two VA pairs aggregated in a four site segment. Further 
examples with a detailed view of the segments, of the 
corresponding relaxed spin structures and of the result¬ 
ing spin currents are reported in Supplementary Note 3 
Supplementary Fig. 5. 

The temperature step of our simulations is 0.4 K 
(0.8 K) for the clean (disordered) case. To better deter¬ 
mine the various phases, at each temperature we con¬ 
struct a histogram over the Monte Carlo history de¬ 
fined on a three-dimensional grid spanned by the or¬ 
der parameter ((/),P(i__i)) (see Ref. [69]). The 
probability for a value (</>, P(i,i), P(i,-i)) of the or¬ 
der parameter is given by the Boltzmann factor ~ 
exp[—P((/), P(i 1 ), P(i _i))/(fcBT)], where F is the free 


energy. Finding the position of the maximum of the 
histogram [which is a point in three-dimensional space 
((/), P(i,i), P(i,_i))] is then equivalent to minimize the free 
energy and identifies the stablest phase. This yields 
sharper transitions than following the thermal average 
which, with our accessible system sizes, often is not large 
enough to resolve closely separated transitions. More de¬ 
tails are given in Supplementary Figs. 9 and 10 and Sup¬ 
plementary Note 5. 
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